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Abstract 

We seek to derive the probability-expressed in terms of the Hilbert-Schmidt (Euclidean or flat) 
metric-that a generic (nine-dimensional) real two-qubit system is separable, by implementing the 
Q^' well-known Peres-Horodecki test on the partial transposes (PTs) of the associated 4x4 density 

' matrices {p). But the full implementation of the test-requiring that the determinant of the PT be 

^ , , nonnegative for separability to hold-appears to be, at least presently, computationally intractable. 

, 

So, we have previously implemented-using the auxiliary concept of a diagonal-entry-parameterized 
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^ ■ separability function (DESF)-the weaker implied test of nonnegativity of the six 2x2 principal 

. minors of the PT. This yielded an exact upper bound on the separability probability of 

I 0.76854. Here, we extend this line of work by requiring that the four 3x3 principal minors of 

' the PT be nonnegative, giving us an improved/reduced upper bound of g « 0.628571. Numerical 

■ simulations-as opposed to exact symbolic calculations-indicate, on the other hand, that the true 



probability is certainly less than i. Our combined analyses lead us to suggest a possible form for 



^ [ the true DESF, yielding a separability probability of || « 0.453125, while the best exact lower 

bound estabhshed so far is 6928 - 22057r ^ 0.0348338. 
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We direct the reader's attention to Fig. [T], which depicts various forms of "diagonal- 
entry-parameterized separabihty functions" (DESF's) jl|, ^-as opposed to "eigenvalue- 
parameterized separability functions (ESFs) jsl, 0, Qj-that we will employ here to obtain 
estimates and simple exact upper bounds on the Hilbert-Schmidt (HS) probability that 
a generic (nine-dimensional) real two-qubit system is separable. Zyczkowski, Horodecki, 
Sanpera and Lewenstein, in a much-cited article Q], have given "philosophical", "practi- 
cal" and "physical" reasons for studying "separability probabilities". We have examined 
the associated problems which arise, using the volume elements of several metrics of in- 
;erest as measures on the quantum states, in various numerical and theoretical studies 
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121. 



The subordinate of the three curves in Fig. [T]-derived using an extensive quasi-Monte 
Carlo (Tezuka-Faure 13, 1^) six-dimensional numerical integration procedure-provides an 



estimate of the true, but so-far not exactly-determined DESF. The dominant of the three 



curves-readily obtainable from results already reported in jl, sec. VII]-has the form 

(3e2? - 1) ^>o 
S^iO = { ■ (1) 

-ie« (e2? - 3) e < 

The intermediate of the three curves, which we first report here, has the same-differing only 
in constants-functional form 



9n^ 



2048 



(7e2« - 27) ^ < 



With each of these three curves we can obtain an associated estimate or upper bound on 
the desired HS separability probability {P^p/reai)- This is accomplished by integrating over 
^ G [—00,00] the product of the corresponding curve with the function (Fig. [2]) (based on 
the jacobian of a coordinate transformation, to be described below) 

64csch^(0(-160 sinh(20 - 25 sinh(40 + 12^(16 cosh(20 + cosh(40 + 18)) 



Psl'/real= I "^lO^lOde. (4) 



that is, 

L 

Sep /real 

Proceeding, thusly, we obtain an upper bound on the HS separability probability of 
~ 0.76854 based on the dominant of the three curves, the titular || ^ 0.628571 
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FIG. 1: Three forms of diagonal-entry-parameterized separability functions (DESFs) 
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FIG. 2: Jacobian ([3]), which when multiplied by a separability function and integrated over ^ G 
[—00,00], yields the associated Hilbert-Schmidt separability probability 

using the intermediate curve, and an estimate of 0.4528427 for the true probability with 
the subordinate, numerically- derived curve. (From our work in [4, eq. (25)], we already 
know that the HS probability of a generic real two-qubit system being absolutely separable- 
that is not entanglable by any unitary transformation-is ^'^'^^^'/2^^^ ~ 0.0348338, which then 
serves as a lower bound on the corresponding HS [absolute plus nonabsolute] separability 
probability itself (cf. Q Q, eq. (29)]).) 

The variable ^ used in the above presentation is the logarithm of the square root of the 
ratio of the product of the 11- and 44-entry of the associated real 4x4 density matrix (p) 
to the product of the 22- and 33-entries, that is 

P11P44 



e = lo^ 
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iia, 



(In our previous studies [l|, l2] , we have employed the alternative variables, v = and 
fi = ^ p^^p33 ; but now switch to the [more symmetric] form ([5]). Importantly, only the "cross- 
product ratio of diagonal entries is needed in our parameterization to test for separability, 
and not the individual entries themselves.) The jacobian (j3]) used in our calculations is 
obtained by the transformation of one of the diagonal entries, say, P33, to ^ and integrating 



the Hilbert-Schmidt (Lebesgue) volume element (of course, P44 = 1 — pn — P22 ^ P33) [la, P- 
13646] 

3/3 

dVnS = (PllP22P33P44)~dpiidp22dp33, /3 = 1 (6) 

over pii and P22 and normalizing the result. (To obtain the corresponding HS volume 
elements for the complex 4x4 density matrices, one must employ-conforming to a pattern 
familiar from random matrix theory-/9 = 2, and [3 = 4 in the quaternionic case (cf. fl?!).) 
The use of the celebrated Peres-Horodecki separability test Q, Q is central to our 



analyses. Ideally, we would be able to require that t 



of p be nonnegative to guarantee separability 



le determinant of the partial transpose 
211 ]. However, this has so far proved to be 
too computationally demanding a (fourth-degree, high-dimensional) task for us to enforce 
(cf. [Ij, eq. (7)]). But, in j^, we did succeed in implementing the weaker implied test that 
all the six 2x2 principal minors of the partial transpose of p be nonnegative, giving us 
the dominant curve in Fig. [H (Actually, only two of the minors differ nontrivially from 
the analogous set of minors of p itself.) To derive the sharper intermediate curve here, we 
extended this approach to the four 3x3 principal minors. Actually, we found that these 
four minors fell into two pairs of identical results. Further, one of the set of results 

f > 

9 (n = .1 (7) 

^3x31^,; ^ 37re-3«(e«v/l^(37e2«+2e4«+2l)+3(27e2«-7)sm-i(e«)) 

1024 < 

could be obtained from the other set by the transformation ^ —>■ — ^. This curve ([7]) and 
its reflection around ^ = are shown in Fig. [3l The intermediate curve ([2]) in Fig. [H first 
reported here, was constructed by joining the sharper segments of these two curves over the 
two half-axes. (A parallel strategy had been pursued with the 2x2 minors.) The comparable 
results to and Fig. [3] for the 2x2 minors investigation [2] are 

e-2«(2 sinh(n + cosh(n) ^ > 
S2.2iO = { (8) 
1 ^ < 
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FIG. 3: The two distinct (red and blue) separability functions obtained from the four 3x3 principal 
minors, the envelope of which defines the intermediate curve in Fig. [TJ The intercept, as well as 
the point of intersection of the two curves, is at ~ 0.867446. 
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FIG. 4: The two distinct (red and blue) separability functions obtained from the six 2x2 principal 
minors, the envelope of which defines the dominant curve in Fig. [1] 



and Fig. HI 

For the intermediate curve in Fig. [1] we have the nontrivial y-axis intercept of 



457r^ 
512 



0.867446 (the intercept for the dominant curve being simply 1), while the estimate of the 
true intercept using the numerically-generated curve is 0.612243, quite close to our previously 



conjectured value of 



1357r^ 
2176 



0.612315 



a. 



In obtaining our several results, we used the " Bloore/correlation" parameterization of 
density matrices 22, |23| and accompanying ranges of integration-generated by the cylin- 



drical algebraic decomposition procedure 2m, l25|, implementing the requirement that p be 
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nonnegative definite-presented in [l], eqs. (3)-(5)]. The computational tractability of utiliz- 
ing the 3x3 principal minors of the partial transpose in this coordinate frame appeared to 
stem from the fact that each of these four quantities only contains three of the six off-diagonal 
variables {zij) employed in the full parameterization (each set of three variables, additionally 
and conveniently, sharing a common row/column subscript). (The nine-dimensional convex 
set of real two-qubit density matrices is parameterized by six off-diagonal-^ = -^^^-and 
three diagonal variables-pjj.) Integrating out the three variables not present in the con- 
straint simply leaves us with a constrained (boolean) integration over the cube [—1,1]^, as 
indicated in eq. (3)]. We appropriately permuted the subscripts in the indicated coordi- 
nate system, so that we could study all four of the minors (thus, finding that they fell into 
two equal sets). Of course, such a simplifying integration strategy is not available for the 
determinant of the partial transpose itself, which contains all the six off-diagonal variables 
(zjj), rather than simply three. 

Each of the constrained integrations we have utilized so far, has used as it constraint 
the nonnegativity of a single 2 x 2 or 3 x 3 principal minor of the partial transpose of p. 
(However, we have been able above to couple and complement multiple results by taking 
the sharper/tighter bounds over the half-axes provided by individual outcomes.) We have, 
to this point, been unable-using either the (Bloore [2j]) density-matrix parameterization 



23|- 



presented in jl| or the interesting partial-correlation parameterization indicated in 
to perform constrained integrations in which two or more 2 x 2 or 3 x 3 minors (and a 
fortiori the determinant) are required to be simultaneously nonnegative. (It, then, remains 
an open question whether or not being able to do so would simply lead to the dominant and 
intermediate curves given already in Fig. [1] and by ([1]) and ([2])). 

We can, however, rather convincingly-but in a somewhat heuristic manner-reduce the 
derived upper bound on the HS separability probability of generic real two-qubit systems 
from II ~ 0.628571 to 0.576219 by using a new curve-having a y-intercept of {^^Y ~ 
0.752462 as a DESF. This curve is obtained by taking the product oi the two curves displayed 
in Fig. [3] (that is, the product of the function ([2]) with its reflection about ^ = 0). A plot of 
the result shows that it is both subordinate to the intermediate curve in Fig. [H as is obvious 
it must be, but also clearly dominates the numerically-generated curve there, which is an 
estimate of the true DESF. (Since each of the two curves in Fig. H] is simply unity over a 
half-axis, a parallel strategy in the 2x2 minors analysis can, of course, yield no nontrivial 
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upper-bound reduction from ~ 0.76854.) 

The "twofold-ratio" theorem of Szarek, Bengtsson and Zyczkowski 



26l|-motivated by 



the numerical results reported in lOl] -allows us to immediately obtain exact upper bounds, 
as well, on the HS separability probability for generic (eight-dimensional) real minimally- 
degenerate real two-qubit systems (boundary states having a single eigenvalue zero). These 
upper bounds would, then, be one-half those applicable to the nondegenerate case-that is, 
^§^2 ^ 0.38427 and § ^ 0.314286. Further, we can, using the results of our numerical 
study, similarly obtain an induced estimate, 0.226421, of the true probability. 

The two sets of derived functions ([T]) and ([2]), based respectively on the 2x2 and 3x3 
minors have the same functional forms, but with differing sets of constants ({1,2,3, 1} vs. 
{9,2048 = 2^^,27,7}). It seems natural, then, to conjecture that the true separability 
'unction-which must be based on the determinant of the partial transpose P, eq. (7)] 



2d . |21| . that is, the single 4x4 minor-will also adhere to the same functional form, but 
with a different set of constants. 

In fact, pursuing this line of thought, we have found that the function 

SlSe-aSf-S+lSeaOvra ^ ^ 

— ^ > 

SconjectureiO = \ 315e« (-18+5e2« (9) 

St6 — ^ < 

fits (Fig. [5]) the numerically-generated subordinate curve in Fig. [1] quite well, yielding an 
HS separability probability of f = 0.453125, and a y-intercept of ^ 0.6167. 



(Then, by the twofold-ratio theorem 26|], the HS separability probability of the minimally- 



degenerate (boundary) states would be ^ = P ^ 0.226563. Also, we have been able to find 
a number of other curves, adhering to this same general structure, fitting the subordinate 
curve in Fig. [T]equally as well, and again yielding || as a separability probability, in addition 



to well-fitting curves yielding somewhat less simple fractions-such as ||| ~^0. 452778, ^ 
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0.453086 and ||| ^ 0.45291.) We are obligated, however, to note that in ^, sec. IX.A] we 
had advanced-based on somewhat different considerations than here-the hypothesis that 
this probability is ^ ~ 0.470588, with an associated DESF equal to 

{135e-3«(-l+3e2e)7r2 
kiT — ^ > 
135e«(-3+e^e).2 ' (^O) 

(However, our best numerical estimate at that point was 0.4538838 pj, sec. V.A.2] [j, sec. 
IX.A] , rather close to our current-study estimate of 0.4528427. By computing standard errors 
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FIG. 5: The difference between the numericaUy-generated subordinate function in Fig. [5] and a 
suggested possibly true separabiUty function ([9]), giving a separabihty probabiUty of || 0.453125 
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FIG. 6: The difference between the numerically-generated subordinate function in Fig. [5] and the 
previously conjectured true separability function (fTO]) . giving a separability probability of ^ fti 
0.470588 

of the mean, we can establish a [~ 95%] confidence range-four standard deviations wide-for 
this estimate of (0.451634, 0.454051)-that does contain || ^ 0.453125. A comparable plot 
(Fig. [6]) to Fig. [5] shows (ITO!) to provide a considerably poorer fit.) 

One might further speculate-in line with random matrix theory and our previous analyses 
j2|-that the DESF for the generic (15-dimensional) complex two-qubit systems is propor- 
tional to the square of ([9]). If the constant of proportionality were simply taken to equal unity, 
the associated HS separability probability would be f^guTeroel = ^''^23?.n^^" ~ 0.252864, 
rather close to the value ^ ~ 0.242424 conjectured, for a number of reasons, in ^, sec. 
IX.B]. 
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